Stripes and the t-J model 
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We investigate the two-dimensional t-J model at a hole doping of x = | and J/t = 0.35 with exact 
diagonalization. The low-energy states are uniform (not striped). We find numerous excited states 
with charge density wave structures, which may be interpreted as striped phases. Some of these are 
consistent with neutron scattering data on the cuprates and nickelates. 
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In the search for an understanding of the cuprate su- 
perconductors it is desirable to find a model which cap- 
tures many of the essential aspects of the environment 
experienced by the electrons in these materials. Because 
these materials are born out of antiferromagnetic insula- 
tors by doping, it is somewhat urgent to decide if a simple 
model which begins from the strong electron correlation 
limit, such as the so-called t-J model, can explain some 
features which result from the electronic degrees of free- 
dom in the cuprates jjj . Even though the progress made 
in trying to solve the t-J model may be characterized as 
slow, it gives some features which are present in these 
materials. For example, some important aspects of the 
calculated single-hole spectrum [|| are in agreement with 
the results of the photo-emission data In addition, 
the model gives rise to a two-hole bound state with 
the d x 2_ y 2 symmetry which is the believed symmetry of 
the superconducting state in these materials. 

Emery and Kivelson suggested that the cuprates 
are near an electronic phase separation instability which 
is prevented by the long-range part of the Coulomb in- 
teraction. In the phase-separated state, the holes cluster 
together, leaving the rest of the system in an antiferro- 
magnet state with no holes. Phase separation in the t-J 
model has been studied by a number of techniques which 
seem to be giving conflicting conclusions [|o]-|f0|. For ex- 
ample, Hellberg and Manousakis ||] using a stochastic 
projection method, an extension of the Green's function 
Monte Carlo (GFMC) for lattice fermions, find that the 
t- J model has a region of phase separation at all interac- 
tion strengths. Other techniques fail to reach this con- 
clusion. Most of these studies use small size systems 0, 
high temperature series expansions or approximate 
methods §,0. In a very recent calculation Calandra et 
al. [[l0| using the GFMC approach within the fixed node 
approximation find that the phase boundary for phase 
separation is far from that determined by the high tem- 
perature series expansions [|| and much closer to that 
obtained by Hellberg and Manousakis [|| except in the 
delicate region with small hole dopings and J/t < 0.4. 
By using a uniform Fermi-liquid type nodal structure 
one disregards the possibility of a non-uniform ground 
state in which one component of the mixture (the an- 



tiferromagnetic phase) has no fcrmion degrees of free- 
dom. In addition, in the delicate region of small J/t 
and low doping, Shraiman and Siggia |fli"| and Bonin- 
segni and Manousakis showed spin-back-flow effects 
become very important resulting in the interesting struc- 
ture of the hole "polaron". These effects are known to 
change the nodal structure of the wave function in a cru- 
cial way in strongly correlated quantum fluids. Therefore 
fixed-node GFMC may be inadequate in this region. 

The phase separation in the t-J model cannot be real- 
ized in the physical system due to the Coulomb interac- 
tion m|b|. Instead such a tendency for phase separa- 
tion can be satisfied locally by forming stripes or other 
charge density wave (CDW) structures without a large 
Coulomb cost. The coupling of electrons to lattice dis- 
tortions may also encourage the formation of stripes. 

Experimentally, stripe modulations were first observed 
in a doped nickelate analogue of the cuprates Jl4) . 
La 2 Ni04 may be doped with holes by adding oxygen or 
by substituting strontium for lanthanum. The modu- 
lation seen with neutron scattering in the doped com- 
pounds is consistent with the holes forming diagonal do- 
main walls separating antiferromagnetic regions of spins. 
Stripes with a variety of widths and hole densities along 
the stripes have been observed. 

There is strong evidence for stripe modulations in the 
cuprates as well In Lai e-aNdo^Sr^CuC^, super- 

conductivity is suppressed at a filling of i = g, and 
neutron scattering studies reveal vertical domain walls 
of holes and spins. In these stripes, a hole density of 
Ph = 1/2 per lattice spacing is observed. 

Recently White and Scalapino (WS) found static ver- 
tical stripe order similar to that of the cuprates in the 
two-dimensional t-J model using a density matrix renor- 
malization group technique p6[ . These results are sur- 
prising due to the fact that the t-J model ignores the 
long-range part of the Coulomb interaction and couples 
to no lattice distortions. One sees no physical reason for 
such a simplified model to have a ground state with a 
periodic array of interfaces. 

In this paper we study the two-dimensional t-J model 
at a hole doping of x = g and at J/t = 0.35 with ex- 
act diagonalization. The low-energy states are uniform 
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(not striped). We find a variety of CDW excited states 
which may be interpreted as striped phases. Some of 
these are in excellent agreement with neutron scattering 
data on the cuprates and the nickelates. However, with- 
out adding additional terms to the t-J model, such as a 
long-range Coulomb interaction or a coupling to the lat- 
tice, the striped states are only realized as excited states. 

The t-J Hamiltonian is written in the subspace with 
no doubly occupied sites as 



mrij 



(1) 



Here (ij) enumerates neighboring sites on a square lat- 
tice, cj ff creates an electron of spin a on site i, m = 
Sct c L c *o"' an d is the spin-i operator. Throughout 
this paper, we take t = 1 and J = 0.35. 

To achieve a doping of a; = |, all calculations were 
carried out on periodic 16-site clusters with two holes. 
Periodic clusters may be characterized by their primitive 
translation vectors ai and &2- There are a large number 
of possible 16-site clusters on the two-dimensional square 
lattice. Each cluster can only support striped phases that 
are commensurate with the periodicity of that particular 
cluster [ pT[ . Clusters which have one particularly short 
translation vector are quasi-one-dimensional and behave 
like chains or ladders. We only consider clusters in which 
both translation vectors have a Manhattan length of at 
least Im — \a x \ + \ aV \ > 4. There are seven such clusters, 
shown in Table |[ These clusters represent all possible 
quasi-two-dimensional 16-site clusters. Later in the pa- 
per, we examine six eigenstates in detail. We label these 
states by letters, (a) through (f), shown next to their 
corresponding clusters in Table [ij. 

TABLE I. Translation vectors of the seven possible 
quasi-two-dimensional 16-site clusters. The letters label 
eigenstates examined later in the paper. 



Number 



ai 



a.- 



States 



(0,4) 
(0,4) 
(0,4) 
(3,2) 
(3,1) 
(2,2) 
(2,2) 



(4,0) 
(4,1) 
(4,2) 
(2,-4) 
(L-5) 
(4,-4) 
(3,-5) 



(e),(f) 



(d) 
(a),(b) 
(c) 



We studied all eigenstates with energy per site £ < 
—0.634 for each of the seven periodic clusters. For 
each cluster, we use all possible combinations of phases 
9 = 0, 7T along each translation vector. Several of these 
states are CDWs, which are necessarily degenerate. The 
density in a CDW state is characterized by an amplitude 
A, a wave- vector k, and an arbitrary phase 4>. The phase 
<t> depends on the particular linear combination of degen- 
erate states taken. Thus the hole density on site is 
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FIG. 1. Amplitude of the CDW in every state with energy 
per site £ < —0.634. The six labeled states correspond to 
states examined in detail later in the paper. The low-energy 
states have no CDWs. 



given by 



n h (Ti) = n h + Acos(k • 



Ti + ■ 



(2) 



where the average hole density is fih = g • 

The amplitude A of the CDW in every low-energy state 
are plotted in Fig. [I] as a function of energy per site. The 
lowest energy states are uniform and have CDW ampli- 
tude A = 0. For energies above £ > —0.645 some CDW 
states are stabilized. The maximum CDW amplitude of 
these states increases with increasing energy. 

We examine in detail the six labeled states in Fig. [I]. 
State (a) is the lowest energy state and is uniform. States 
(b) through (f) have increasing CDW amplitude and in- 
creasing energy. Each of these states has the largest 
CDW amplitude of all states at or below its energy. 

We examine the charge order of the six states in Fig. ||. 
The CDW states, (b) through (f), are degenerate, so tak- 
ing a different linear combination of the eigenstates will 
move the CDW. In particular, in state (b) the maximum 
charge order occurs on a site, while in states (c), (d), (e), 
and (f), the maximum charge order occurs between two 
sites, or on a bond. All of the states have site-centered 
and bond-centered CDWs with different linear combina- 
tions. 

States (b) and (c) have diagonal stripes with two dif- 
ferent hole densities along the stripe. State (b) has hole 
density pb = 1/2 per (1,1) step, while state (c) has p c = 1. 
These states are similar to experimental results on the 
nickelates [Q and to the mean-field calculations of Zaa- 
nen and Littlewood fjsf ]. 

States (e) and (f) exhibit vertical stripes with 1/2 hole 
per (0,1) step consistent with experimental results on 
Lai.6- a; Ndo.4Sr a ;Cu04 |15j and with the calculations of 
WS (l§. States (e) and (f) have the largest CDW am- 
plitudes that we found. 

Interestingly, the CDWs of the vertical stripe states, 
(e) and (f), are remarkably similar to the density profile 
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FIG. 2. Hole structures of the six eigenstates we chose to 
examine in detail. The radius of each circle is proportional 
to the hole density on the given site. Additionally, the cir- 
cles are shaded according to the relative hole density in each 
state: Black circles show the maximal hole density in that 
state while white circles show the minimum. State (a), which 
has the minimum energy, has a uniform hole density. States 
(b) through (f) have increasing energy and show increasing 
CDW amplitude. In all states the average spin moment on 
each site is zero. 
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FIG. 3. Comparison of the hole density rij, as a function of 
^-coordinate for a 4 x 4 cluster with open boundary conditions 
in the ^-direction and states (e) and (f). All three systems 
have exhibit a stripe in the y-direction. The state with open 
boundary conditions is the non-degenerate ground state. 



obtained from the ground state of a 4 x 4 cluster with 
open boundary conditions in the a;-direction and periodic 
boundary conditions in the the y-direction, the boundary 
conditions used by WS |l6| ]. The hole density as a func- 
tion of the ^-coordinate for these three states is shown in 
Fig. 0. The states (e) and (f ) are excited states of the pe- 
riodic cluster #2 in Table |. The 4x4 cluster with open 
boundary conditions in the ^-direction can be generated 
from the periodic cluster #2 by cutting all bonds along 
one column. Clearly, the excited states (e) and (f) essen- 
tially have an extra node along one column. The open 
boundary conditions in one direction causes the ground 
state of this cluster to be very similar to excited states 
of the periodic cluster. The open boundary conditions 
select the striped state to be the ground state. 

None of the eigenstates has a spontaneous spin density 
wave amplitude, but the spin correlations are affected by 
the CDW. One way to show the spin correlations is to 
apply small magnetic fields along the boundary of the 
simulation cell, as done by WS jl6). This is shown for 
the diagonal and vertical stripe states (b) and (f) in Fig. 
[|. The stripes are pinned so the sites with the fields have 
maximum electron density with the appropriate average 
polarization. In both the diagonal and vertical cases, the 
antifcrromagnetic order in neighboring stripes is shifted 
by 7T, as in the nickelates and cuprates |l4|]l5| ]. 

To conclude, we found stripes in the t- J model at a dop- 
ing of x = |, but only as excited states. The ground state 
of the model for J/t = 0.35 is uniform. The energy cost 
per site to form diagonal stripes similar to those found 
in the nickelates is at least Ad > 0.016t, and for vertical 
stripes similar to those in the cuprates the energy cost 
is A v > 0.0252;. At this doping and interaction strength 
studies on larger systems have found that the model is 
nearly phase separated ||. The stripes seen experimen- 
tally could be the result of phase separation frustrated 
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FIG. 4. Hole and spin structures of states (b) and (f) gen- 
erated by applying a small staggered magnetic field of h = 0.1 
to sites at the boundary of the simulation cell to pin the spins. 
The radius of each circle is proportional to the hole density 
on the given site, and the length of the arrows is proportional 
to {3^). These are not eigenstates of the t-J Hamiltonian 
without the pinning fields, but the spin correlations pictured 
are similar to the correlations in the eigenstates (b) and (f). 
State (e) has a spin and charge structure very similar to state 
(f) and is not shown. 



by the Coulomb repulsion J13[ and/or the coupling of 
the electrons to lattice distortions. However, the stripe 
states are not ground states of the simple t-J model, as 
claimed in Ref. fll6[ . 

Our findings are in agreement with the recent work of 
Pryadko, Kivelson and Hone Jl9| who studied the interac- 
tion between localized holes in a weakly doped quantum 
antiferromagnet. They find that stripes are unstable due 
to an attractive interaction between such domain walls. 
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at Wright Patterson Air Force Base. 
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